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ABSTRACT 

We present z band photometry of three consecutive transits of the exoplanet 
TrES-1, with an accuracy of 0.15% and a cadence of 40 seconds. We improve 
upon estimates of the system parameters, finding in particular that the planetary 
radius is 1.081 ± 0.029 R Jnp and the stellar radius is 0.811 ± 0.020 R Q . The 
uncertainties include both the statistical error and the systematic error arising 
from the uncertainty in the stellar mass. The transit times are determined to 
within about 15 seconds, and allow us to refine the estimate of the mean orbital 
period: P = 3.0300737 ± 0.0000026 days. We find no evidence for star spots or 
other irregularities that have been previously reported. 

Subject headings: planetary systems — stars: individual (TrES-1) — techniques: 
photometric 

1. Introduction 

The passage of a planet in front of its parent star is an occasion for celebration and 
for intensive observations. Transits provide a wealth of information about an exoplanetary 
system, even when there is little hope of ever resolving the system with adaptive optics, 
coronagraphy, or interferometry (see, e.g., the recent review by Charbonneau et al. 2006). 
The dimming events can be recorded photometrically (as first done by Charbonneau et 
al. 2000 and Henry et al. 2000) and used to determine the planetary and stellar radii. Short- 
term anomalies in the timing of the transits may betray the presence of moons or other 
planets (Holman & Murray 2005, Agol et al. 2005), and long-term variations in the light- 
curve shape should result from orbital precession (Miralda-Escude 2002), although neither 
effect has yet been observed. A few stellar absorption lines have been observed to deepen 
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due to absorption by the planetary atmosphere (Charbonneau et al. 2002, Vidal-Madjar et 
al. 2003). There are also anomalous Doppler shifts due to the Rossiter-McLaughlin effect, 
which reveal the angle on the sky between the stellar rotation axis and the orbital axis (Queloz 
et al. 2000; Winn et al. 2005,2006; Wolf et al. 2006). Secondary eclipses — the counterpoints 
to transits, when the planet passes behind the star — have not yet been detected at optical 
wavelengths, where the signal would reveal the planetary albedo. However, they have been 
detected at mid-infrared wavelengths (Charbonneau et al. 2005; Deming et al. 2005, 2006) 
and are beginning to reveal details of exoplanetary thermal emission. 

The aim of the Transit Light Curve (TLC) project is to build a library of high-precision 
transit photometry. Our main scientific goals are to support all of the transit investigations 
mentioned above, by refining estimates of the basic system parameters; and to search for 
secular and short-term variations in transit times and light-curve shapes that would be 
produced by as-yet undiscovered planets or moons. We have previously reported on TLC 
observations of the transiting exoplanets XO-lb (Holman et al. 2006) and OGLE-TR-lllb 
(Winn et al. 2006). 

In this paper, we present TLC results for TrES-lb, whose discovery by Alonso et 
al. (2004) was notable for being the first success among the many ground-based, wide- 
field surveys for transiting planets with bright parent stars. The TrES-1 parent star is 12th 
magnitude with spectral type K0 V. The planet is a "hot Jupiter" with a mass and radius 
of approximately 0.8 Mj up and 1.0 i?j up , and an orbital period that is almost exactly 3 days. 
Calculations by Laughlin et al. (2005) have shown that the measured planetary radius is in 
accordance with models of irradiated hot Jupiters. Charbonneau et al. (2005) have detected 
thermal emission from the planet. Steffen & Agol (2005) have searched for transit timing 
anomalies, with null results. 

This paper is organized as follows. We describe the observations in the next section, 
and the photometric procedure in § 3. In § 4, we describe the techniques we used to estimate 
the physical and orbital parameters, and in § 5 we provide the results, along with a closer 
look at the characteristics of the photometric noise. A brief summary is given in § 6. 

2. Observations 

We observed three transits of TrES-1 (on UT 2006 June 9, 12, and 15), corresponding 
to epochs E = 234, 235, and 236 of the ephemeris given by Alonso et al. (2004): 



T C (E) = 2,453, 186.8060 [HJD] + E x (3.030065 days). 



(1) 
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We used KeplerCam on the 1.2m (48 in) telescope at the Fred L. Whipple Observatory 
(FLWO) on Mt. Hopkins, Arizona. KeplerCam has one 4096 2 Fairchild 486 back-illuminated 
CCD, with a 23'1 x 23'1 field of view. For our observations we used 2x2 binning, which gives 
a scale of 0'/68 per binned pixel, a readout /setup time of 11 s, and a typical readout noise of 
7 e~ per binned pixel. We observed through the SDSS z band filter in order to minimize the 
effect of color-dependent atmospheric extinction on the relative photometry, and to minimize 
the effect of limb-darkening on the transit light curve. The effective band pass was limited 
at the blue end by the filter (which has a transmission of 25% at 825 nm, rising to 97% for 
wavelengths redward of 900 nm) and at the red end by the quantum efficiency of the CCD 
(which drops from 90% at 800 nm to 25% at 1000 nm). 

On each of the three nights, we observed TrES-1 for approximately 6 hr bracketing the 
predicted midpoint of the transit. We acquired 30 second exposures in good focus. We also 
obtained dome flat exposures and zero-second (bias) exposures at the beginning and the end 
of each night. 

On the night of UT June 9, we observed under generally light clouds, through an 
airmass ranging from 1.01 to 1.35. A period of thicker clouds caused about 45 minutes 
of data following egress to be unusable. There were also small temperature fluctuations 
(~0.1 deg min -1 ) that proved to cause systematic errors in the photometry, as described in 
§ 3. In addition, the guider was not functioning as intended, causing the position of TrES-1 
on the detector to drift by about 20 pixels over the course of the observations. The full- width 
at half- maximum (FWHM) of stellar images ranged from 2 to 3 binned pixels (1'.'3 to 2'/0). 

The weather was nearly perfect on the nights of UT June 12th and 15th, with no visible 
clouds. Automatic guiding was functioning, and the pixel position of TrES-1 varied by no 
more than 4 pixels over the course of each night. On the 12th, the airmass ranged from 
1.01 to 1.18, and the FWHM was nearly constant at approximately 2.4 binned pixels (1'.'6). 
On the 15th, the airmass ranged from 1.01 to 1.20, and the FWHM ranged from 2.5 to 
3.5 binned pixels (l'/7 to 2'.'3). 
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3. Data Reduction 

We used standard IRAF0 procedures for the overscan correction, trimming, bias sub- 
traction, and flat-field division. We performed aperture photometry of TrES-1 and 12 nearby 
stars of comparable brightness. Many different aperture sizes were tried in order to find the 
one that produced the minimum noise in the out-of-transit data. For the three nights, the 
optimum aperture radii were 7, 8, and 7 pixels (4'.'7, 5'.'4, and 4f7). We subtracted the under- 
lying contribution from the sky, after estimating its brightness within an annulus centered 
on each star ranging from 30 to 35 pixels in radius. 

For each night's data, the following steps were followed to produce relative photometry 
of TrES-1: (1) The light curve of each comparison star was normalized to have unit median. 
(2) A comparison signal was created by taking the mean of all 12 normalized light curves, 
after rejecting obvious outliers in the comparison star light curves. (3) The light curve of 
TrES-1 was divided by the comparison signal. Although the transit was obvious, there was 
a small residual time gradient in the out-of-transit data, a sign of a residual systematic 
error. (4) To correct for the time gradient and set the out-of-transit flux to unity, the light 
curve was divided by a linear function of time. This function was determined as part of the 
model-fitting procedure and will be described in the next section. 

For the data from UT June 9, an intermediate step between steps (3) and (4) was 
necessary. Those data suffered from an additional systematic error, namely, a few "bumps" 
in the light curve, with amplitude ~0.2% and time scale ~15 minutes. The deviations 
between this light curve and the light curves based on data from the subsequent nights 
proved to be well correlated with both the ambient temperature (especially the temperature 
measured at the mirror cell), and the measured ellipticity of stellar images. Presumably the 
temperature fluctuations were affecting the focus, which in turn caused differential leakage 
of light outside the digital apertures. To correct for this systematic error, we used the data 
from the other two nights to construct a light-curve model (see § 4), and then found the 
residuals between the UT June 9 data and the model. A linear function provided a good fit 
to the correlation between the residuals and ellipticity; we used this relation to correct each 
data point based on the measured ellipticity. We verified that the resulting light curve had 
residuals that were uncorrelated with temperature, FWHM, and ellipticity. On the other 
nights, neither the temperature nor the ellipticity varied as much. The data did not have 
any noticeable correlations with external variables, and no corrections were applied. 



1 The Image Reduction and Analysis Facility (IRAF) is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities for Research in Astronomy, Inc., under 
cooperative agreement with the National Science Foundation. 
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To estimate the relative uncertainties in our photometry, we computed the quadrature 
sum of the errors due to the Poisson noise of the stars (both TrES-1 and the comparison 
stars), the Poisson noise of the sky background, the readout noise, and the scintillation noise 
(as estimated according to the empirical formulas of Young 1967 and Dravins et al. 1998). 
The dominant term is the 0.1% Poisson noise from TrES-1. For the UT June 9 data, an 
additional error term was included in the sum, equal to one-half the size of the ellipticity- 
based correction. 

The final photometry is given in Table 1, and is plotted in Fig. [TJ The quoted uncer- 
tainties have been rescaled by a factor specific to each night, such that xV-^Vdof = 1 when 
fitting each night's data individually (see § 4). The scale factors by which the calculated 
errors for each night needed to be increased were 1.17, 1.13, and 1.26. A composite light 
curve, produced by time-shifting the first and last nights' data appropriately, is shown in 
Fig. [2j The lower panel shows a time-binned version of the composite light curve, in which 
the bin size is 3 minutes. 



4. The Model 

To estimate the planetary, stellar, and orbital parameters, we fitted a parameterized 
model to all of the photometry simultaneously. The model is based on a star and a planet on 
a circular orbit about the center of massEI The star has a mass M* and radius R+, and the 
planet has a mass M p and radius R p . The orbit has a period P and an inclination i relative 
to the sky plane. We define the coordinate system such that 0° < i < 90°. It is often useful 
to refer to the impact parameter b = acosi/R* (where a is the semimajor axis) rather than 
the inclination. 

We allowed each transit to have an independent value of T c , the transit midpoint, rather 
than forcing them to be separated by exact multiples of the orbital period. This is because 
we sought to measure or bound any timing anomalies that may indicate the presence of 
moons or additional planets in the system. Thus, the period P was relevant to the model 
only through the connection between the total mass and the orbital semimajor axis. We 
fixed P = 3.030065 days, the value determined by Alonso et al. (2004). The uncertainty of 
0.000008 days was negligible for this purpose (although we were able to use the resulting 
values of T c to refine the period estimate, as described in § 5). 



2 We assume a circular orbit because, in the absence of any evidence for additional bodies in the system, 
it is expected that tidal forces between the star and the planet have had sufficient time to circularize the 
orbit (see, e.g., Rasio et al. 1996, Trilling et al. 2000, Dobbs-Dixon et al. 2004). 



3 1.00 



> 0.99 



K 0.98 




0.75 



0.80 



HJD - 



0.85 

2453895 



0.90 



0.95 



3 1.00 



> 0.99 



£ 0.98 



UT 2006 June 12 





C_> 
I 

o 



0.80 



0.85 0.90 
HJD - 2453898 



0.95 



3 1.00 



> 0.99 



5 0.98 



UT 2006 June 15 




0.80 



0.85 



0.90 
HJD - 2453901 



ill,,!,! 




0.95 



1.00 



Fig. 1. — Relative z band photometry of TrES-1. The best-fitting model is shown as a solid 
line. The residuals (observed — calculated) and the rescaled 1 a error bars are also shown. 
The residuals have zero mean but are offset by a constant flux to appear beneath each light 
curve, for clarity. From top to bottom, the root-mean-squared (RMS) residuals are 0.16%, 
0.14%, and 0.15%. 
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Fig. 2. — Composite transit light curve of TrES-1. Top. Data from UT June 9, 12, and 15 
are shown with red, green, and blue points, respectively. The best-fitting model is shown as 
a solid line. The residuals (observed — calculated) and the 1 a error bars are shown beneath 
each light curve. The mean time between points is 15 seconds, and the RMS residual is 
0.15%. Bottom. Same, after binning in time with 3 minute bins. The RMS residual is 
0.046%. 
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The stellar mass cannot be determined from transit photometry alone. Furthermore, the 
values of R* and R p that result from fitting the photometry depend on the choice of stellar 
mass§| Our approach was to fix M* at a value that has been estimated from an analysis of 
the stellar spectrum (its luminosity class, effective temperature, surface gravity, etc.) and 

1/3 1/3 

theoretical isochrones, and then use the scaling relations R p oc M* and R+ oc MJ to 
estimate the systematic error due to the uncertainty in M*. We adopted the value M* = 
0.89 ± 0.05 M Q based on the analysis by Sozzetti et al. (2004), noting that this value is also 
in agreement with independent work by Laughlin et al. (2005) and Santos et al. (2006). The 
planetary mass M p is nearly irrelevant to the model (except for its minuscule effect on the 
semimajor axis), but for completeness we use the value M p = 0.76 M Jup found by Sozzetti 
et al. (2004). 

To calculate the relative flux as a function of the projected separation of the planet and 
the star, we assumed the limb darkening law to be quadratic, 

^ = l_ Ul (l_ A1 )_ U2 (l_ M )2 j (2) 

where I is the intensity, and [i is the cosine of the angle between the line of sight and the 
normal to the stellar surface. We employed the analytic formulas of Mandel & Agol (2002) to 
compute the integral of the intensity over the portion of the stellar disk hidden by the planet. 
There is a quandary over how to treat the limb darkening parameters when fitting transit 
photometry. Most previous investigators in this situation have either fixed the parameters 
at the values calculated with stellar atmosphere models (see, e.g., Jha et al. 2000, Moutou 
et al. 2004, Laughlin et al. 2005), or allowed them to vary freely (e.g., Brown et al. 2001, 
Deeg et al. 2001, Winn et al. 2005; see also Knutson et al. 2006, who tried both approaches 
on the same data). However, the first option is problematic because the stellar atmosphere 
models are not perfect, and it may be a mistake to trust them so completely. On the other 
hand, the second option probably allows too much freedom in the limb darkening function, 
and hence too much uncertainty in the fitted parameter values. 

Our approach was to allow u\ and «2 to be free parameters, but with a mild a priori 
constraint based on the stellar atmosphere models of Claret (2004). In particular, we added 
a penalty term to x 2 ( see below) to encourage the limb-to-center intensity ratio Iq/Ii to 
agree within 20% of the value calculated by Claret (2004) based on the ATLAS quadratic 



3 For a fixed P, an increase in the stellar mass causes a to increase in proportion to MJ . This in turn 
causes the orbital velocity of the planet to be increased and the time scale of the transit to be decreased by 

1/3 

the same factor. This can be compensated for by increasing both i?* and R p in proportion to MJ , with 
no other observable effect. 
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fit (u\ = 0.2851, U2 = 0.2849). We chose the figure of 20% based on a survey of the 
recent literature in which limb darkening calculations are compared to interferometric and 
microlensing observations. There seems to be agreement to within 5-20% in the limb-to- 
center intensity ratio for various stellar types (see, e.g., Abe et al. 2003, Fields et al. 2003, 
Ohishi et al. 2004, Aufdenberg et al. 2005, Bigot et al. 2006). We also investigated the effect 
of either tightening or dropping this a priori constraint, as discussed below. 

As noted in the previous section, the light curves exhibited a time gradient in the 
out-of-transit data, probably due to differential extinction between the target star and the 
comparison stars, or some other systematic error. For this reason, each transit was described 
with two additional parameters: the out-of-transit flux / oot and a time gradient a. 

In total, there were 14 adjustable parameters describing 1149 photometric data points. 
The parameters were R+, R p , and i; the limb-darkening parameters Ui and m 2 ; and the values 
of T c , f oot and a for each of the 3 transits. In practice, it was preferable to use the parameters 
v\ = 2ui + ui and v 2 = u\ — 2w 2 rather than u\ and w 2 . This is because v i and t> 2 have 
nearly uncorrelated uncertainties. 

We optimized the parameters by using the AMOEBA algorithm (Press et al. 1992, p. 
408) to minimize the error statistic 



where /j(obs) is the flux observed at time j, o~j is the corresponding uncertainty, and /j(calc) 
is the calculated value. The second term is the a priori constraint on the limb darkening 
function. As noted in § 3, the uncertainties o~j were the calculated uncertainties (based on 
Poisson noise, scintillation noise, etc.), after multiplication by a factor specific to each night, 
such that y 2 /Ndoy = 1 when each night's data was fitted individually. (Our intention was 
not to test the model, but rather to determine the appropriate weights for the data points.) 
The optimized model is plotted as a solid line in Figs. [1] and [2j The differences between the 
observed fluxes and the calculated fluxes are also shown beneath each light curve. 

The statistical uncertainties in the fitted parameters were estimated using a Markov 
Chain Monte Carlo (MCMC) simulation (for a brief introduction, consult appendix A of 
Tegmark et al. 2004). In this method, a stochastic process is used to create a sequence 
of points in parameter space that approximates the joint probability distribution of the 
parameter values, given the data values. The sequence, or "chain," is generated from an 
initial point by iterating a "jump function" which can take various forms. Our jump function 
was the addition of a Gaussian random number to each parameter value. If the new point has 
a lower \ 2 than the previous point, the jump is executed; if not, the jump is only executed 




+ 



(Vjo) - 0.43 
0.086 



(3) 
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with probability exp(— A% 2 /2). We set the perturbation sizes such that ~25% of jumps 
are executed. We created 20 independent chains, each with 500,000 points, starting from 
random initial positions. The first 100,000 points were not used, to minimize the effect of 
the initial condition. The correlation lengths were ^100 for the highly covariant parameters 



R p , i?*, and b, and ~200 for the other parameters. The iGelman fc Rubin! (119921 ) R statistic 



was within 0.2% of unity for each parameter, a sign of good mixing and convergence. 

Table 2 gives the results for each parameter, along with some useful quantities derived 
from the parameters. Among the latter are the times between first and last contact (tjv — 
t{) and between first and second contact (in — ti), which are useful for planning future 
photometric measurements, and the quantity (R p /a) 2 , which controls the amount of starlight 
that is reflected from the planet. 

Fig-Elshows the probability distributions for the especially interesting parameters R+, R p 
and b, along with some of the 2-d probability distributions for some correlated parameters. 
Although the distributions shown in Fig. [3] are somewhat asymmetric about the median, 
Table 2 reports only the median p med and a single number a p characterizing the statistical 
error. The value of a p is the average of |p mc d — Phil and |p me d — Pio|> where p\ Q and p^ are the 
lower and upper 68% confidence limits. Two-sided formal uncertainties are of little practical 
significance for R p and R* because those parameters are subject to a systematic error that 
is comparable to the statistical error, namely, the covariance with the stellar mass. The size 
of this systematic error is also given in Table 2. The results for b and i are best regarded as 
one-sided limits, since central transits are allowed with reasonable probability. 



5. Results 

Our results for the planetary and stellar radii are R p = 1.081 ± 0.029 i?j up and R+ = 
0.811±0.020 Rq. The quoted error in each parameter is the quadrature sum of the statistical 
error (based on the 68% confidence limits found through the MCMC simulations) and the 
systematic error (based on the 1 a uncertainty in the stellar mass quoted by Sozzetti et 
al. 2004). The statistical and systematic errors are comparable in size. The ratio of radii 
is independent of the stellar mass and is known with greater precision: R p /R± = 0.13686 ± 
0.00082. 

Previous investigators have not had photometry of sufficiently high cadence and signal- 
to-noise ratio to solve for both R* and R p simultaneously. Instead, they have used estimates 
for R* (as well as M*) based on an analysis of the stellar spectrum and theoretical isochrones. 
Sozzetti et al. (2004) found R± = 0.83 ± 0.05 R Q , and Laughlin et al. (2005) found R± = 
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Fig. 3. — Top row. Probability distributions for the stellar radius R±, planetary radius R p , 
and impact parameter b = a cos i/R*, based on the MCMC simulations. The arrows mark the 
values of the parameters that minimize x 2 . A solid line marks the median of each distribution, 
and the dashed lines mark the 68% confidence limits. For b, the dashed line marks the 95% 
confidence upper limit. Middle and bottom rows. Joint probability distributions of 
those parameters with the strongest correlations. The contours are isoprobability contours 
enclosing 68% and 95% of the points in the Markov chains. The dotted lines indicate the 
values of the limb darkening parameters calculated by Claret (2004) using an ATLAS model 
of a star with T cff = 5250 K, logg = 4.5, [Fe/H] = 0.0, and & = 1.0 km s _1 . 
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0.83 ±0.03 R Q n Our photometry has allowed us to estimate R* independently, subject only 
to the fairly weak covariance with M*. Our result of R* = 0.811 ± 0.020 R Q confirms the 
previous estimates and improves upon their precision. 

The planetary radius of R p = 1.081±0.029 -Rj up is also in agreement with (and more pre- 
cise than) the previous estimates of 1. 041^05 -^Jup (Sozzetti et al. 2004) and 1.08 ±0.05 i?j up 
(Laughlin et al. 2005). TrES-1 is an important case study for theoretical models of the struc- 
ture of hot Jupiters, in part because of the comparison to the intensively studied exoplanet 
HD 209458. The radius of TrES-1 is smaller by ~20%, despite a similar mass and a similar 
degree of stellar insolation. Calculations by Laughlin et al. (2005) and Baraffe et al. (2005) 
have shown that the properties of TrES-1 are generally in line with theoretical expectations, 
and that HD 209458 is anomalous. This provides the important clue that any mechanism to 
"inflate" HD 209458 cannot apply to all hot Jupiters of the same mass. Showman & Guil- 
lot's (2002) proposal, that atmospheric circulation patterns deposit a significant fraction of 
the stellar insolation into the planetary interior, provides no apparent reason why TrES-1 and 
HD 209458 should be so different, unless perhaps the circulation patterns depend strongly 
on surface temperature (~1060 K for TrES-1 and 1130 K for HD 209458). Bodenheimer 
et al. (2001) proposed that HD 209458 is heated internally by eccentricity tides. Winn & 
Holman (2005) proposed that HD 209458 is trapped in a Cassini state with nonzero obliq- 
uity, causing ongoing heating through obliquity tides. This is naturally a rare occurrence, 
the implication being that most hot Jupiters (including TrES-1) have a negligible obliquity. 
Adding to the puzzle, two newly discovered exoplanets also appear to be anomalously large: 
HAT-P-1 (Bakos et al. 2006) and WASP-1 (Collier Cameron et al. 2006, Charbonneau et 
al. 2006). 

We also confirm that the transit is fairly central, i.e., the trajectory of the planet comes 
close to the center of the stellar disk. With 95% confidence, b < 0.28 and i > 88?4. All 
other things being equal, central transits have the maximum possible duration and allow for 
the greatest number of in-transit measurements (and therefore the greatest signal-to-noise 
ratio). This makes them favorable for follow-up photometry. However, they are unfavorable 
for certain applications, such as the measurement of orbital precession through changes in the 
impact parameter (Miralda-Escude 2002), or of spin-orbit alignment through the Rossiter- 
McLaughlin effect (Ohta et al. 2005, Gaudi & Winn 2006). Specifically, a central transit 
provides little leverage on the determination of the angle A between the projected orbital 
axis and the projected stellar spin axis, because in that case A is strongly covariant with the 
stellar rotation rate. 



4 Laughlin et al. (2005) also made use of the existing photometry showing that the transits are nearly 
central (b « 0) to place a lower limit of 0.80 Rq on the stellar radius. 
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The fitted values of the limb darkening parameters are U\ = 0.28 ± 0.06 and U2 = 
0.21 ± 0.12, with a strong covariance as noted previously. In comparison, the values cal- 
culated by Claret (2004) are U\ = Ui = 0.29. These are in good agreement, which is not 
surprising because of the a priori constraint on the limb darkening function. If this constraint 
is dropped, then the optimized limb darkening parameters are u\ = 0.36 and U2 = 0.01 (im- 
plying a limb-to-center intensity ratio of 63% instead of the theoretical value of 42%), and 
the radius estimates become R+ = 0.815 i?j up and R p = 1.096 -Rj up - Thus, by dropping the a 
priori constraint, the favored values of the radii are driven to larger values, but by less than 
the 1 a statistical error determined previously. Conversely, if the limb darkening parameters 
are held fixed at the Claret (2004) values, the new radius estimates are R+ = 0.810 R@ and 
R p = 1.072 -Rj U p, again within the 1 a statistical error of our quoted results. Limb darkening 
should therefore be considered as another source of systematic error, but one that is no larger 
than (and is probably smaller than) the systematic error due to the covariance with stellar 
mass. 

The statistical errors in our transit times range from 12 to 17 seconds. We used the 
ephemeris of Alonso et al. (2004) to compute "observed minus calculated" (O— C) residuals 
for the newly measured transit times, along with 14 of the transit timea^ and the time of 
secondary eclipse reported by Charbonneau et al. (2005). The residuals are shown in Fig. [H 
We fitted a straight line to the transit times as a function of epoch number E to derive a 
refined ephemeris 

T C (E) = T c (0) + E x P, (4) 

finding T c (0) = 2,453,186.80603(28) [HJD] and P = 3.0300737(13) days, where the numbers 
in parentheses are the uncertainties in the last 2 digits. It must be noted that a line is not 
a statistically acceptable fit to the timing data; the value of \ 2 is 23.5 with 12 degrees of 
freedom (xV^dof = 2). For this reason, we increased all of the errors by a factor of \pl 
before determining the uncertainties in T c (0) and P. The poor fit could be a signal that 
the period is not constant, as assumed, or it could be a signal of systematic errors in the 
measured times. To be conservative, we recommend basing future predictions on the period 
3.0300737 days but with an uncertainty of 0.0000026 days, i.e., twice as large as the formal 
uncertainty. Interestingly, our 3 new measurements are only marginally consistent with a 
uniform period at the 1 a level. The transits occur progressively later than expected. Of 
course, with only 3 measurements this could be a coincidence, but future measurements are 
motivated by the possible detection of a satellite or additional planet through a periodicity 



5 We excluded the 2 timing measurements that were based on observations of only part of the transit, 
i.e., for which data were missing either prior to ingress or after the egress. This is because pre- ingress and 
post-egress data are very important for assessing and correcting systematic errors in the photometry. 
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Fig. 4. — Transit and secondary-eclipse timing residuals for TrES-1, from Charbonneau et 
al. (2005) and this work. The calculated times (using the ephemeris of Alonso et al. 2004) 
have been subtracted from the observed times. The best-fitting line is plotted, representing 
the updated ephemeris given in Eq. (TjJ and the text that follows it. 

in the residuals. 

Previous investigators have reported on interesting photometric anomalies in TrES-1 
light curves. The review paper by Charbonneau et al. (2006) includes an exquisite light 
curve obtained with the Hubble Space Telescope in which there is a positive flux excursion 
relative to the model of ~0.2%, lasting for ^30 minutes (Brown et al., in preparation). The 
most natural interpretation is that the planet occulted a star spot or group of star spots, 
a phenomenon anticipated by Silva (2003). We find no evidence for star spots in our data, 
and can set bounds on spot properties under certain assumptions. For a circular spot much 
smaller than the planet, the anomaly in the transit light curve caused by the spot occultation 
is approximately 



where -R spo t is the spot radius, J spot is the intensity of the spot, and is the intensity of the 
unspotted stellar photosphere in the vicinity of the spot. The intensity ratio can be related 
to the temperature difference via the blackbody formula, as in Eq. (1) of Silva (2003). The 
duration of sunspot-eclipse totality is ~2_R p /w orb , or 20 minutes, if the spot lies exactly on 
the transit chord. A more likely impact parameter of 0.5 would lead to a duration of about 
15 minutes. By averaging our data into 15 minute time bins, we find random residuals with 
a standard deviation of o = 2 x 10~ 4 . Thus we can set an approximate 2a upper limit of 
/ < 4 x 10~ 4 on the product of the areal ratio and intensity contrast of small, circular star 
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spots. This upper bound is too weak to rule out Sun-like starspots (/ ^ 7 x 10~ 6 )Jj although 
the spots on later-type stars like TrES-1 are expected to be more prominent than sunspots 
(see, e.g., Alekseev 2006). We can rule out flux anomalies as large as the one seen by Brown 
et al., but there is no real contradiction because the observations were taken over a year 
apart. 

In addition, an enterprising group of amateur astronomers created a composite light 
curve of TrES-1 transits based on over 10,000 individual brightness measurements (Price 
et al. 2006). They reported statistically significant evidence for a "brightening episode" by 
0.5% during egress. We find no evidence for such episodes in any of our light curves, which 
have 0.15% accuracy per 40-second interval. 

To characterize our noise statistics more generally, the left panel of Fig. shows a 
histogram of the flux residuals (observed — calculated) from all three nights of data. The 
distribution of residuals is approximately Gaussian with a standard deviation of 0.15%. 
The residuals are not noticeably correlated with the pixel position of TrES-1, the ambient 
temperature, or the shape parameters of stellar images. We also time-averaged the composite 
light curve using time bins of various sizes, and calculated the standard deviation of the 
residuals as a function of the time bin size t. The results are shown in the right panel of 
Fig. [5j The noise is reduced as 1/ \ft over two orders of magnitude, from the unbinned mean 
cadence of 15 seconds up to 30 minutes, at which point we do not have enough data to 
usefully average further. 

6. Summary 

Through observations of three consecutive transits, we have significantly improved upon 
the estimates of the system parameters for TrES-1. Our results are in agreement with previ- 
ous results, but are more precise. We thus confirm the previous conclusions that the planetary 
radius is in general agreement with theoretical predictions for irradiated hot Jupiters. We 
have also made an independent determination of R+ from the photometric signal that con- 
firms the previous estimates that were based on spectral typing and theoretical isochrones. 
We have measured the 3 transit times with an accuracy of about 15 seconds, leading to a 
refined estimate of the mean period, and providing an anchor for continuing searches for 
timing anomalies. 



6 Here, we have taken (R spo t/R*) ~ 10 5 , T spot ~ 4240 K and T+ ss 6050 K based on the typical size and 
temperature contrast of a solar sunspot umbra, from Steinegger et al. (1990). 
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Fig. 5. — Noise properties of the flux residuals. To compute residuals, the model-predicted 
relative flux has been subtracted from each point of the composite transit light curve shown in 
Fig. [2] Left. The distribution of unbinned residuals. The dotted line is a Gaussian function 
with a standard deviation of 0.15%. Right. The standard deviation of the residuals, as 
a function of the size of the time- averaging bin size. The solid line represents the l/yt 
dependence that is expected in the absence of systematic errors. 

The noise in our photometry is nearly Gaussian and averages down with the expected 
1/ \/t dependence all the way to an averaging time of 30 minutes. We achieved similar results 
for the exoplanet XO-1 (Holman et al. 2006), although in that case we measured only two 
transits with the FLWO 1.2m telescope and Keplercam. We find this apparent randomness 
of the residuals to be very encouraging, as it raises the possibility of detecting signals at the 
~10 -4 or even smaller levels, such as those produced by smaller transiting planets, moons, 
and reflected light. These projects are generally thought to be the exclusive province of 
space-based satellite photometry, but it may be possible to achieve them through extensive 
and repeated observations with a relatively small ground-based telescope. 

We thank Scott Gaudi for very helpful comments about our analysis and the manuscript. 
We thank Perry Berlind and Mike Calkins for help with the observing; Max Tegmark for help 
with the MCMC simulations; and Dave Charbonneau, Ben Lane, and Dimitar Sasselov, for 
useful discussions. A.R. thanks the MIT UROP office for research funding. Keplercam was 
developed with partial support from the Kepler Mission under NASA Cooperative Agreement 
NCC2-1390 (PI D. Latham). 
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Table 1. Photometry of TrES-1 



HJD 


Relative flux 


Uncertainty 


2453901.96003 


1.00029 


0.00154 


2453901.96049 


0.99700 


0.00154 


2453901.96094 


1.00009 


0.00154 


2453901.96139 


0.99985 


0.00154 


2453901.96185 


0.99998 


0.00154 


2453901.96230 


0.99894 


0.00154 


2453901.96275 


1.00039 


0.00154 


2453901.96320 


1.00139 


0.00154 



Note. - The time stamps represent the 
Heliocentric Julian Date at the time of mid- 
exposure. The uncertainty estimates are based 
on the procedures described in § 2. We intend 
for this Table to appear in entirety in the elec- 
tronic version of the journal. A portion is shown 
here to illustrate its format. The data are also 
available from the authors upon request. 



- 21 - 



Table 2. System Parameters of TrES-1 



Parameter 


Value 


Uncertainty 


R*/Rq 


0.811 


0.020 a 


Rpl Rjup 


1.081 


0.029 a 


W0.89 Mq)- 1 ^ 


0.811 


0.012 


(i? P /i?jup)(M./0.89 M )-V3 


1.081 


0.021 


Rp/R* 


0.13686 


0.00082 


R*/a 


0.0957 


0.0014 


(Rp/a) 2 


0.000171 


0.000006 


i [deg] 


> 88.4 


(95% conf.) b 
(95% conf.) b 


6 


< 0.28 


irv - k [hr] 


2.497 


0.012 


*n — ii [min] 


18.51 


0.63 


Ul 


0.284 


0.061 




0.21 


0.12 


2ui + u 2 


0.777 


0.046 


U\ — 2u 2 


-0.13 


0.29 


T c (234) [HJD] 


2453895.84297 


0.00018 


T c (235) [HJD] 


2453898.87341 


0.00014 


T c (236) [HJD] 


2453901.90372 


0.00019 



Note. — The values in Column 2 are the medians p mo d of the 
MCMC distributions, which are based on the assumption M* = 
0.89 M Q and an a priori constraint on the limb darkening function 
(see the text). Except where noted, the value in Column 3 is 
the average of |p mod - Pi Q | and \p mcd - p hi \, where p lo and p hi 
are the lower and upper 68% confidence limits. (The cumulative 
probability for p < p\ Q is 16%, and the cumulative probability for 
p > phi is 16%.) 

a These quantities are affected by a systematic error associated 
with the stellar mass. The quoted uncertainties are the quadra- 
ture sum of the statistical error and the systematic error (esti- 
mated under the assumption M+ — 0.89 ± 0.05 M ). 

b Given the shapes of the probability distributions for b and i, 
these results are best regarded as one-sided limits. 



